Spin transport and tunable Gilbert damping in a single-molecule magnet junction 
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We study time-dependent electronic and spin transport through an electronic level connected 
to two leads and coupled with a single-molecule magnet via exchange interaction. The molecular 
spin is treated as a classical variable and precesses around an external magnetic field. We derive 
expressions for charge and spin currents by means of the Keldysh non-equilibrium Green's functions 
technique in linear order with respect to the time-dependent magnetic field created by this precession. 
The coupling between the electronic spins and the magnetization dynamics of the molecule creates 
inelastic tunneling processes which contribute to the spin currents. The inelastic spin currents, in 
turn, generate a spin-transfer torque acting on the molecular spin. This back-action includes a 
contribution to the Gilbert damping and a modification of the precession frequency. The Gilbert 
damping coefficient can be controlled by the bias and gate voltages or via the external magnetic 
field and has a non-monotonic dependence on the tunneling rates. 

PACS numbers: 73.23.-b, 75.76.+J, 85.65.+h, 85.75.-d 



I. INTRODUCTION 

Single-molecule magnets (SMMs) are quantum mag- 
nets, i.e. mesoscopic quantum objects with a per- 
manent magnetization. They are typically formed by 
paramagnetic ions stabilized by surrounding organic 
ligands. 1 SMMs show both classical properties such 
as magnetization hysteresis^ and quantum properties 
such as spin-tunnelingP coherence^ and quantum phase 
interferen ce^^. They have recently been in the center 
of intereslPSH in view of their possible applications as 
information storage^ and processing devices P 

Currently, a goal in the field of nanophysics is to 
control and manipulate individual quantum systems, in 
particular individual spins EHUD Some theoretical works 
have investigated electronic transport through a molecu- 
lar magnet contacted to leads.ESHSDln this case, the trans- 
port properties are modified due to the exchange inter- 
action between the itinerant electrons and the SMMj 2 ^ 
making it possible to read out the spin state of the 
molecule using transport currents. Conversely, the spin 
dynamics and hence the state of an SMM can also be 
controlled by transport currents. Efficient control of the 
molecule's spin state can be achieved by coupling to fer- 
romagnetic contacts as well. 21 

Experiments have addressed the electronic transport 
pro pertie s through magnetic molecules such as Mni2 and 
p eg j22]23] w hj c h have been intensively studied as they 
are promising candidates for memory devices!^ Various 
phenomena such as large conductance gaps, 2 ^ switch- 
ing behavior J2§1 negative differential conductance, depen- 
dence of the transport on magnetic fields and Coulomb 
blockade have been experimentally observed! 22 * 23 * 27 * ^ 
Experimental techniques, inc luding fo r instance scanning 
tunneling microscopy (STM), 22 23 29 3 3 break junctions^ 
and three-terminal devices, ^ 2 * 23 * 27 * have been employed 
to measure electronic transport through an SMM. 
Scanning tunneling spectroscopy and STM experiments 



show that quantum properties of SMMs are preserved 
when deposited on substrates!^ The Kondo effect in 
SMMs with magnetic anisotropy has been i nvestigated 
both theoretically^ anc j experimentallyP^H it has been 

suggested^ and experimentally verified^ that a spin- 
polarized tip can be used to control the magnetic state 
of a single Mn atom. 

In some limits, the large spin S of an SMM can be 
treated as a classical magnetic moment. In that case, 
the spin dynamics is described by the Landau-Lifshitz- 
Gilbert (LLC) equation that incorporates effects of ex- 
ternal magnetic fields as w ell as torques originating from 
damping phenomena! 3 ' 3S ' In tunnel junctions with mag- 
netic particles, LLG equat ions have been derived using 
perturbative couplings^HSI and the nonequilibrium Born- 
Oppenheimer approximation.^ Current-induced magne- 
tization switching is driven by a generated spin-transfer 
torque (STT^M^H^Jas a back-action effect of the electronic 
spin transport on the magnetic particle ft 6 ! 45 ^ ^ A spin- 
polarized STMTS' has been used to experimentally study 
STTs in relation to a molecular magnetizationPS This 
experimental achievement opens new possibilities for 
data storage technology and applications using current- 
induced STTs. 

The goal of this paper is to study the interplay be- 
tween electronic spin currents and the spin dynamics of 
an SMM. We focus on the spin transport properties of 
a tunnel junction through which transport occurs via a 
single electronic energy level in the presence of an SMM. 
The electronic level may belong to a neighboring quan- 
tum dot (QD) or it may be an orbital related to the 
SMM itself. The electronic level and the molecular spin 
are coupled via exchange interaction, allowing for inter- 
action between the spins of the itinerant electrons tunnel- 
ing through the electronic level and the spin dynamics of 
the SMM. We use a semiclassical approach in which the 
magnetization of the SMM is treated as a classical spin 
whose dynamics is controlled by an external magnetic 
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field, while for the electronic spin and charge transport 
we use instead a quantum description. The magnetic 
field is assumed to be constant, leading to a precessional 
motion of the spin around the magnetic field axis. The 
electronic level is subjected both to the effects of the 
molecular spin and the external magnetic field, generat- 
ing a Zeeman split of the level. The spin precession makes 
additional channels available for transport, which leads 
to the possibility of precession-assisted inelastic tunnel- 
ing. During a tunnel event, spin-angular momentum may 
be transferred between the inelastic spin currents and the 
molecular spin, leading to an STT that may be used to 
manipulate the spin of the SMM. These torques include 
the so-called Gilbert damping, which is a phenomenolog- 
ically introduced damping term of the LLG equation, 38 
and a term corresponding to a modification of the pre- 
cession frequency. We show that STTs and hence the 
SMM's spin dynamics can be controlled by the external 
magnetic field, the bias voltage across the junction and 
the gate voltage acting on the electronic level. 

The paper is organized as follows: We introduce our 
model and formalism based on Keldysh non-equilibrium 
Green's functions technique^^ in Sec. [ilj where we 
derive expressions for the charge and spin currents in lin- 
ear order with respect to a time-dependent magnetic field 
and analyze the spin transport properties at zero temper- 
ature. In Sec. Ill we replace the general magnetic field of 
Sec. [II] by an SMM whose spin precesses in an external 
constant magnetic field, calculate the STT components 
related to the Gilbert damping and the modification of 
the precession frequency and analyze the effects of the 
external magnetic field as well as the bias and gate volt- 
ages on the spin dynamics. Conclusions are given in Sec. 
IV. 



II. CURRENT RESPONSE TO A TIME 
DEPENDENT MAGNETIC FIELD 

A. Model and Formalism 

For the sake of clarity, we start by considering a junc- 
tion consisting of a non-interacting single-level QD cou- 
pled with two normal, metallic leads in the presence of an 
external, time-dependent magnetic field (see Fig. [TJ . The 
leads are assumed to be noninteracting and unaffected 
by the external field. The total Hamiltonian describing 
the junction is given by H(t) = Hl,r + H T + H D (t). 
The Hamiltonian of the free electrons in the leads reads 
H L ,R = J2k,a,^e{LM} e kaic\a^k^, where £ denotes the 
left (L) or right (R) lead, whereas the tunnel coupling 
between the QD and the leads can be written as Ht = 
Efe.a^ei.-R^e^Lf^ + V ka^Uk^}- The tunnel ma- 
trix element is given by Vfc CT £. The operators $is a p{ckoi) 

and a^ida) are the creation (annihilation) operators of 
the electrons in the leads and the QD, respectively. The 
subscript a = J [, J, denotes the spin-up or spin-down state 




FIG. 1: (Color online) A quantum dot with a single electronic 
level, eo, coupled to two metallic leads with chemical poten- 
tials /j, l and fiR in the presence of an external time-dependent 
magnetic field, B(t). The spin-transport properties of the 
junction are determined by the bias voltage eV = fix, — [ir, 
the position of the level eo, the tunnel rates I\ and Tr and 
the magnetic field. 

of the electrons. The electronic level, e , of the QD is 
influenced by an external magnetic field, B(t). consisting 
of a constant part, B c , and a time-dependent part, B'(t). 
The Hamiltonian of the QD describing the interaction be- 
tween the electronic spin s and the magnetic field is then 
given by Hr)(t) — Hfj + H'(t), where the constant and 
time-dependent parts are = 'Yl llT £od\ r d a + g/j,BsB c 
and H'(t) — gussB' (t). The proportionality factor g is 
the gyromagnetic ratio of the electron and /is is the Bohr 
magneton. 

The average charge and spin currents from the left lead 
to the electronic level are given by 

I L u{t) = q v {j t N L »] =q v ^([H,N Lv ]), (1) 

where N Lv = 2fe, CT , CT / ^iO^W' 2 *^ is the charge and 
spin occupation number operator of the left contact. The 
index v — corresponds to the charge current, while 
v = x,y, z indicates the different components of the spin- 
polarized current. The current coefficients q v are then 
go = — e and q v ^o = h/2. In addition, it is useful to 
define the vector o v = (1, a) , where I is the identity 
operator and a consists of the Pauli operators with ma- 
trix elements (c?) CTO -'. Using the Keldysh nonequilibrium 
Green's functi ons technique, the currents can then be ob- 
tained atPES 

M*) = - J dt'Tr{a v [G r (t,t')±<(t',t) (2) 

+ G<(t,t')± a L (t',t)}}, 

where G r,a,< are the retarded, advanced and lesser 
Green's functions of the electrons in the QD with the ma- 
trix elements G^",(t,t') = =fi6{±t T t')U^(*). 
and G< a ,(t,t') = i(S a ,(f)d a {t)), while Ej ,< (t,t') are 
self-energies from the coupling between the QD and 
the left lead. Their matrix elements are given by 
{^<{t,t')) a , a = E.VW^'^MO^W The 
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Green's functions gj,'^ (t,t') are the retarded, advanced 
and lesser Green's functions of the free electrons in the 
left lead. The retarded Green's functions G r of the elec- 
trons in the QD, in the presence of the constant mag- 
netic field, B c . are found using the equation of motion 
technique,^! while the lesser Green's functions, Gq , are 
obtained from the Keldysh equation Gq — GqY, < Gq, 
where multiplication implies internal time integrations.^ 
The time-dependent part of the magnetic field can be 

expressed as B'(t) = E w (^ e_i "' + K eiut )> where 
is a complex amplitude. This magnetic field acts as a 
time-dependent perturbation that can be expressed as 
it'{t) = J2^(Huje^ luJt +Hle' ut ), where H u is an operator 
in the electronic spin space and its matrix representaton 
in the basis of eigenstates of s z is given by 



9HB 



B^JX "F i-B^y 



(3) 



Applying Dyson's expansion, analytic continuation rules 
and the Keldysh equation,^ one obtains a first-order ap- 
proximation of the Green's functions describing the elec- 
trons in the QD that can be written as 



(4) 



The expression for the currents in this linear approxima- 
tion is given by 



2g„ 



ReTr{a4G r £<+G<X a L 



(5) 



H G 1, L + G H G L L + G H G 2, L J). 



Eq. (|5| is then Fourier transformed in the wide- 
band limit, in which the level width function, T(e) = 
— 2Im{£ r (e)}, is constant, Re{S r (e)} = and one can 
hence write the retarded self-energy originating from the 
dot-lead coupling as E r ' a (e) = =pir/2. From this trans- 
formation, one obtains 



hu(t) = It +^[/ fo (o;)e- iwt + Il v (u) t 



(G) 



Using units in which h = 1, the dc part of the 
currents^ if^ and the time-independent complex com- 
ponents Il v {u}) are given by 

It=1u 1 ^^[/ £ (e)-/«(e)]TtIm{^(e)} (7) 



and 



de T L T R 



{[f L (e) - f R (e)}: 



(8) 



2vr T 

Tr{a,[G r (e - u)H^(e) + 2i hn{G r {e)}H u G a {e + lj)}} 



Z=L,R 



±[f L (e + w) - / e (e)] Tr[a,GS(e)^ w Gg(e + w)]}. 



In the above expressions, /f(e) = [e^~ tl ^/ kBT + is 
the Fermi distribution of the electrons in lead £, where 
kg is the Boltzmann constant. The retarded Green's 
function Gg(e) is given by Gg(e) = [e — eo — 57(e) — 
(l/2) gMs( ?B c ]- 1 P 

The linear response of the spin current with respect 
to the applied time-dependent magnetic field can be ex- 
pressed in terms of complex spin current susceptibilities 
defined as 



dB ui ' 



J = x,y,z. 



(9) 



The complex components Il v {oj) are conversely given by 
Ibvi 1 ^) — Xvji^Buj. By taking into account that 

dH u /dB u j = (l/2)gpiBcrj and using Eq. (JsJ) , the current 
susceptibilities can be written as 



= -<*fi / ^{[/^-/^x (10) 

Tr{a v [G r (e - u)<r 3 G r (e) + 2ilm{G r (e)}a j G%(e + u)}} 



The components obey Xuji -00 ) = Xuj( u )- I n other 
words, they satisfy the Kramers-Kronig relations^! that 
can be written in a compact form as 



x£» = -p 



00 xij(0 



(ii) 



with P denoting the principal value. 

For any k = x,y,z such that j ^ k and j,k 7^ 
i, where z indicates the direction of the constant part 
of the magnetic field B c = B c Si, the complex current 
susceptibilities satisfy the relations 



and 



(12) 
(13) 



in addition to Eq. (111. The other non-zero compo- 
nents are Xoi( a; ) ano - Xui 10 )- ^ n the absence of a constant 
magnetic field, the only non-vanishing components obey 



X 



L (w) = x^(^) = xLH- 



Finally, the average value of the electronic spin in the 

QD reads s(t) = (%)) = (1/2) J2 at7 , ^<4(t)&'(*)> = 

— (i/2) X)o-ct' ®<r"'G'^i a {t,t) and the complex spin suscep- 
tibilities are defined as 



9si(w) 



(14) 



They represent the linear responses of the electronic spin 
components to the applied time-dependent magnetic field 
and satisfy the relations similar to Eqs. (Ill, (12) and 



( 13 1 given above. 
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B. Analysis of the spin and current responses 

We start by analyzing the transport properties of the 
junction at zero temperature in response to the exter- 
nal time-dependent magnetic field, B(t). The constant 
component of the magnetic field, B c , generates a Zeeman 
split of the QD level eo resulting in the levels e^, where 
e-^j, = e ±giiBB c /2 in this section. The time-dependent 
periodic component of the magnetic field, B'(t), then cre- 
ates additional states, i.e. sidebands, at energies e+ ± w 
and e^ ± w (see Fig. |2]). These Zeeman levels and side- 
bands contribute to the elastic transport properties of the 
junction when their energies lie inside the bias-voltage 
window of eV = fJ-L — fi>R- 

However, energy levels outside the bias- voltage window 
may also contribute to the electronic transport due to in- 
elastic tunnel processes generated by the time-dependent 
magnetic field. In these inelastic processes, an electron 
transmitted from the left lead to the QD can change its 
energy by u> and either tunnel back to the left lead or 
out into the right lead. If this perturbation is small, as is 
assumed in this paper where we consider first-order cor- 
rections, the transport properties are still dominated by 
the elastic, energy-conserving tunnel processes that are 
associated with the Zeeman levels. 

The energy levels of the QD determine transport prop- 
erties such as the spin current susceptibilities and the 
spin susceptibilities, which are shown in Fig. [3] The 
imaginary and real parts of the susceptibilities are plot- 
ted as functions of the frequency lo in Figs. [3] (a) and (c). 
In this case, the position of the unperturbed level eo is 
symmetric with respect to the Fermi surfaces of the leads 
and a peak or step in the spin current and spin suscepti- 
bilities appears at a value of w, for which an energy level 
is aligned with one of the lead Fermi surfaces. In Figs. 
[3] (b) and (d) , the susceptibilities are instead plotted as 
functions of the bias voltage, eV. Here, each peak or 
step in the susceptibilities corresponds to a change in the 
number of available transport channels. The bias volt- 
age is applied in such a way that the energy of the Fermi 
surface of the right lead is fixed at hr — while the en- 
ergy of the left lead's Fermi surface is varied according 
to fiL — eV. 



III. SPIN-TRANSFER TORQUES AND 
MOLECULAR SPIN DYNAMICS 

A. Model with a precessing molecular spin 

Now, we apply the formalism of the previous section 
to the case of resonant tunneling through a QD in the 
presence of a constant external magnetic field and an 
SMM; see Fig. [| (a). An SMM with a spin S lives in 
a (25 + l)-dimensional Hilbert space. We assume that 
the spin S of the SMM is large and neglecting the quan- 
tum fluctuations, one can treat it as a classical vector 



_ 
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FIG. 2: (Color online) Sketch of the electronic energy levels of 
the QD in the presence of a time-dependent magnetic field. In 
a static magnetic field, the electronic level eo (solid black line) 
splits into the Zeeman levels e^ (solid red and blue lines). If 
the magnetic field in addition to the static component also in- 
cludes a time-dependent part with a characteristic frequency 
u, additional levels appear at energies e-f ±w (dotted red lines) 
and 64. ± uj (dotted blue lines). Hence, there are six channels 
available for transport. 



whose end point moves on a sphere of radius S. In the 
presence of a constant magnetic field, B c = B c e z , the 
molecular spin precesses around the field axis according 
to S(t) = S± cos(u>i,t)e x + S± sm(u)L,t)e v + S z e z , where 
Sj_ is the projection of S onto the xy plane, ujl = 9^bB c 
is the Larmor precession frequency and S z is the projec- 
tion of the spin on the z axis; see Fig. [4] (b). The spins 
of the electrons in the electronic level are coupled to the 
spin of the SMM via the exchange interaction J. The 
contribution of the external magnetic field and the pre- 
cessional motion of the SMM's spin create an effective 
time-dependent magnetic field acting on the electronic 
level. 

The Hamiltonian of the system is now given by H (t) = 
Hl,r + Ht + Hd (t) + Hs, where the Hamiltonians Hl,r 
and Ht are the same as in Sec. [n] The Hamiltonian 
Hs — gusSB represents the interaction of the molecu- 
lar spin, S, with the magnetic field B c and consequently 
does not affect the electronic transport through the junc- 
tion, but instead contributes to the spin dynamics of the 
SMM. The Hamiltonian of the QD in this case is given by 
H D (t) = H C D +H'(t). Here, H C D - £ CT e dld a +g^sB^ 
is the Hamiltonian of the electrons in the QD in the pres- 
ence of the constant part of the effective magnetic field, 
given by B£ s = [B c + jn^S z \ e z . The second term of the 

QD Hamiltonian, H'(t) — g/iBsB' eS (t), represents the in- 
teraction between the electronic spins of the QD, s, and 
the time-dependent part of the effective magnetic field, 
given by B' cS (t) = j^[cos(uj L t)e x + sm(uj L t)e y ] . The 
time-dependent effective magnetic field can be rewritten 
as B' eB (t) = B LUL e- lu ' Lt +5* £ e iW£ *, where B UJL consists 
of the complex amplitudes B U>LX = JS±/2gfj,B, B ULV — 
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FIG. 3: (Color online) (a) Frequency and (b) bias-voltage dependence of the spin current susceptibilities, (c) Frequency and 
(d) bias- voltage dependence of the spin susceptibilities. In (a) and (c), the chemical potential of the left lead is — 2eo, while 
in (b) and (d), the frequency is set to uj — 0.16eo. All plots are obtained at zero temperature with B c = B c e z , and the other 
parameters set to (ir = 0, ef = 1.48eo, £\ = 0.52eo, T = 0.02eo and Tl ~^r = O.Oleo. 



iJS±/2g/iB and B U>LZ = 0. The time-dependent pertur- 
bation can then be expressed as H'(t) — H UL e~' lu}Lt + 
H} 0L e luJLt , where H UL is an operator that can be written, 
using Eq. (3) and the above expressions for B ULi , as 

JS± ( 1 



Hu L - 







(15) 



The time-dependent part of the effective magnetic field 
creates inelastic tunnel processes that contribute to the 
currents. The in-plane components of the inelastic spin 
current fulfill 



JS ± 



(16) 



2g/j, 



B 



N + iXxyi u >L)], 



where B c is replaced by B^a . The z component vanishes 
to lowest order in iJ'(i). 54 Therefore, the inelastic spin 
current has a polarization that precesses in the xy plane. 
The inelastic spin currents, in turn, exert an STT^S^ 
on the molecular spin given by 



T(t) = I L (t) + I R (t), 



(17) 



thus contributing to the dynamics of the molecular spin 
through 

S(t) = g[isB c x S(t) + T(t). (18) 
Using expressions drjl, (|8| and (15 1, the torque of Eq. 



( 17 1 can be calculated in terms of the Green's functions 



G r (e) and Gg(e) as 

JS± f de sr-^ T^T\ 



2tt ^ r 



\fi(e + u L )-f x {e)]x (19) 



Im{(^)4.t G o,tt( e ) G SrW( e + WL ) £ 



-iuJLt 



}■■ 



with A = L,R. Here ((Ti)if, Ggtt( e )' anc ^ ^o^( e ) are 
matrix elements of <jj, G r a (e) and Gg(e) with respect to 
the basis of eigenstates of s z . This STT can be rewritten 
in terms of the SMM's spin vector as 

f(t) = x S(t) + (3§(t) + 7 S(t). (20) 

The first term in this back-action gives a contribu- 
tion to the Gilbert damping, characterized by the Gilbert 
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(a) 

Ml 




FIG. 4: (a) (Color online) Resonant tunneling in the presence 
of an SMM and an external, constant magnetic field. The 
electronic level of Fig. [I] is now coupled with the spin of 
an SMM via exchange interaction with the coupling constant 
J. The dynamics of the SMM's spin, S, is controlled by the 
external magnetic field, B c , that also affects the electronic 
level, (b) Precessional motion of the SMM's spin in a constant 
magnetic field, B c , applied along the z axis. 



damping coefficient a. The second term acts as an effec- 
tive constant magnetic field and changes the precession 
frequency of the spin, S, with the corresponding coeffi- 
cient f3. The third term cancels the z component of the 
Gilbert damping term, thus restricting the STT to the 
xy plane. The coefficient of the third term, 7, is related 
to a by 7/a = lulS\/ SS Z . Expressing the coefficients 
a and (3 in terms of the current susceptibilities Xxx^l) 
and Xxv^l) results in 



JS Z 



^[Rc{xLK)l - lm{xi v (uL)}}, (21) 



P = ~ EtMxLK)} + Re{ X | y K)}]. (22) 



By inserting the explicit expressions for G r ^{e) and 
Gq (e + ujl), one obtains 



Is J ^£ r * r A[/«(e + ^)-A(e)]x (23) 

uj\ - (e-f - ej 2 
[(L ) 2 + {e _ €tm r )2 + {e _ ei+UJL)2 y 



Try , J^r c r A [/ c ( e + c L )-/ A ( e )]x (24) 

L J £A 



(£) 2 + (e-e t )(e- ei + WL ) 



[(!)2 + (e- et )2] 



) 2 + {(--e i + u: L f 



where = e ±gfj lB B^ s /2 = e ± (u>l + JS z )/2 are the 
energies of the Zeeman levels in this section. In the small 
precession frequency regime, u>l <C fcsT, 7 — > and in 
the limit of S z /S — >• 1 the expression for the coefficient 
a is in agreement with Ref. 1161 



- e; - 
Mb = 



FIG. 5: (Color online) Sketch of the electronic energy levels 
of the QD in the presence of a molecular spin precessing with 
the frequency ojl around an external, constant magnetic field. 
The corresponding Zeeman levels are e-^j,. The precessional 
motion of the molecular spin results in absorption (emission) 
of energy corresponding to a spin flip from spin up (down) to 
spin down (up). Hence, there are only four channels available 
for transport. 



B. Analysis of the spin-transfer torque 

In the case of resonant tunneling in the presence of 
a molecular spin precessing in a constant external mag- 
netic field, one also needs to take the exchange of spin- 
angular momentum between the molecular spin and the 
electronic spins into account in addition to the effects 
of the external magnetic field. Due to the precessional 
motion of the molecular spin, an electron in the QD ab- 
sorbing (emitting) an energy lul also undergoes a spin 
flip from spin up (down) to spin down (up), as is indi- 
cated by the arrows in Fig. [5] As a result, the levels at 
energies e-^j, =F ujl are forbidden and hence do not con- 
tribute to the transport processes. Consequently, there 
are only four transport channels, which are located at 
energies e-^j, ± wl- Also in this case, there are elastic 
and inelastic tunnel processes. Some of the possible in- 
elastic tunnel processes are shown in Fig. [6j These re- 
strictions on the inelastic tunnel processes are also vis- 
ible in Fig. [3] (b), which identically corresponds to the 
case of the presence of a precessing molecular spin with 
uil = 0.16eo and JS Z = 0.8eo- Namely, from Eq. ( |16| , 
which is equivalent to Hb{Ilx(^l)} = L y {uj — 
7S MMxLK)} - Im{ X £> L )}] and lm{I Lx ("L)} = 



-Be{I Ly (uj L )} = 



JS, 



Im{x. 



2 SM l""Uii\-L/j 1 Re{xx y (u L )}], 
and from the symmetries of the susceptibilities displayed 
in Fig. [3](b), it follows that there are no spin currents at 
eV = 6f4 T^l- 

As was mentioned, the spin currents generate an STT 
acting on the molecular spin. A necessary condition 
for the existence of an STT, and hence finite values of 
the coefficients a and f3 in Eqs. (23) and (24), is that 
T L (t) £ -Tr{£) (see Eq. (17l). This condition is met 
by the spin currents generated e.g. by the inelastic tun- 
nel processes shown in Figs. [6] (b) and (c) . These tunnel 
processes occur when an electron can tunnel into the QD, 
undergo a spin flip and then tunnel off the QD into ei- 
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(a) 



(d) 



(b) 



(c) 



(f) 



-liR 



-flR 



FIG. 6: (Color online) Sketch of the inelastic spin-tunneling 
processes in the QD in the presence of the precessing molec- 
ular spin in the field B c = B c e z for different positions of 
the energy levels with respect to the chemical potentials of 
the leads, fiL and {ir. Only transitions between levels with 
the same color (blue or red) are allowed. Different colored 
curved arrows (magenta, brown or green) represent different 
processes. 



ther lead. From these tunnel processes, it is implied that 
the Gilbert damping coefficient a and the coefficient f3 
can be controlled by the applied bias or gate voltage as 
well as by the external magnetic field. If a pair of QD 
energy levels, coupled via spin-flip processes lie within 
the bias-voltage window, the spin currents instead fulfill 
Ji(i) = In.(t), leading to a vanishing STT; see Fig. [(3](d). 
In Figs. [6] (e) and (f) the position of the energy levels of 
the QD are symmetric with respect to the Fermi levels 
of the leads, /i^ and /ir. When the QD level with energy 
e^ + lol is aligned with fi^, this simultaneously corre- 
sponds to the energy level e^ — ljl being aligned with 
Hr; see Fig. [6](f). As a result, a spin-down electron can 
now tunnel from the left lead into the level e^ +ujl, while 
a spin-up electron in the level e^ — ul can tunnel into the 
right lead. These additional processes enhance the STT 
compared to that of the case [6] (e) . 

The two spin-torque coefficients a and /3 exhibit a non- 
monotonic dependence on the tunneling rates, T, as can 
be seen in Figs. [7j H and [9] For T -» 0, it is obvious 
that a, ft — > 0. In the weak coupling limit T <C lol- the 
coefficients a and (3 are finite if the Fermi surface energy 
of the lead £, fulfills either of the conditions 



or e^ < /i£ < + wi 



(25) 
(26) 



in such a way that each condition is satisfied by the 
Fermi energy of maximum one lead. These conditions 
are relaxed for larger tunnel couplings as a consequence 
of the broadening of the QD energy levels, which is also 
responsible for the initial enhancement of a and /3 with 
increasing T. Notice however that a and (3 are eventually 
suppressed for r 3> ojl, when the QD energy levels are 



significantly broadened and overlap so that spin-flip pro- 
cesses are equally probable in each direction and there 
is no net effect on the molecular spin. Physically, this 
suppression of the STT can be understood by noticing 
that for r > Wi a current-carrying electron perceives the 
molecular spin as almost static due to its slow precession 
compared to the electronic tunneling rates and hence the 
exchange of angular momenta is reduced. With increas- 
ing tunneling rates, the coefficient /3 becomes negative 

before it drops to zero causing the torque f3S to oppose 
the rotational motion of the spin S. 

In Fig. [7J the Gilbert damping coefficient a and the 
coefficient (3 are plotted as functions of the applied bias 
voltage at zero temperature. We analyze the case of the 
smallest value of T (red lines), assuming that ojl > 0. 
For small eV, all QD energy levels lie outside the bias- 
voltage window and there is no spin transport; see Fig. 
6 (a). Hence a, (3 —> 0. At eV = e; - Wi the tunnel 
processes in Fig. 6 (b) come into play, leading to a finite 
STT and the coefficient a increases while the coefficient 
(3 has a local minimum. In the voltage region specified by 
Eq. (25) for the coefficient a approaches a constant 
value while the coefficient /? increases. Increasing the 
bias voltage to eV = e^ the tunnel processes in Fig. 6 (c) 
occur, leading to a decrease of a and a local maximum 
of (3. For < eV < €f, the coefficients a, (3 — > 0; see Fig 
6 (d). In the voltage region specified by Eq. (26) for /it/,, 
a approaches the same constant value mentioned above 
while /3 decreases between a local maximum at eV = 
e-f and a local minimum at eV = + u>l, which have 
the same values as previously mentioned extrema. With 
further increase of eV . all QD energy levels lie within the 
bias- voltage window and the STT consequently vanishes. 

Fig. [8] shows the spin-torque coefficients a and /3 as 
functions of the position of the electronic level eo- STTs 
acting on the molecular spin occur if the electronic level 



eo is positioned in such a way that the inequalities ( 25 1 
and ( |26[ ) may be satisfied by some values of eV, es and 
ljl. Again, we analyze the case of the smallest value of 
r (red curve). For the particular choice of parameters in 
Fig. [8] there are four regions for which the inequalities 
p5| and (26 1 are satisfied. Within these regions, a ap- 
proaches a constant value while f3 has a local maximum 
as well as a local minimum. These local extrema occur 
when one of the Fermi surfaces is aligned with one of the 
energy levels of the QD. For other values of e , both a 
and f3 vanish. 

The coefficients a and f3 are plotted as functions of the 
precession frequency lol in Fig. [9j Here, eo = ey/2 and 
therefore the positions of the energy levels of the QD are 
symmetric with respect to the Fermi levels of the leads, 
Hl and fiR. Once more, we focus first on the case of the 
smallest value of T (indicated by the red curve). The 
energies of all four levels of the QD depend on oj^, i.e 
B c . For ojl > 0, when the magnitude of the external 
magnetic field is large enough, the tunnel processes in 
Fig. [6] (f) take place due to the above-mentioned sym- 
metries. These tunnel processes lead to a finite STT, 



0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 2.5 

eV[e ] eV[e ] 



FIG. 7: (Color online) (a) Gilbert damping coefficient a and (b) coefficient j3 as functions of the applied bias voltage eV = 
Mi — Mfli with jj,R — 0, for different tunneling rates F, at zero temperature. Other parameters are Tl = = T/2, e-f = 1.48eo, 
ei = 0.52eo, S = 100, J = O.Oleo, JS Z = 0.8eo and o;l = 0.16eo- In the case of the smallest value of F (red lines), a approaches 
a constant value when lies within the energy range specified by (25) and (26) . The coefficient /3 has one local minimum and 
one local maximum for the same energy range. 




eoteV] e [eV] 

FIG. 8: (Color online) (a) Gilbert damping coefficient a and (b) coefficient j3 as functions of the position of the electronic 
level eo for different tunneling rates F at zero temperature. The applied bias voltage is eV = \il — Hr, with fj,R — 0. Other 
parameters are = Tr — T/2, e t - e = 0.24eV, S = 100, J = 0.005eF, JS Z = OAeV and tv L = 0.08eV. In the case of the 
samllest value of V (red lines), there are four regions in which the Gilbert damping and the change of the precession frequency 
occur. In each of these regions eo satisfies the inequalities (25) and (26), and a approaches a constant value, while j3 has one 
local maximum and one local minimum. 



a maximum for the Gilbert damping coefficient a and 
a negative minimum value for the f3 coefficient. As ujl 
increases, the inequalities of Eqs. (25 1 and (26 1 are sat- 
isfied and the tunnel processes shown in Fig. [dj (e) may 
occur. Hence, there is a contribution to the STTs, but 
as is shown in Eq. ( 23 1 , the Gilbert damping decreases 



with increasing precession frequency. At larger values of 
u>l, resulting in e-f = /zl, the Gilbert damping coefficient 
drops to zero, while the coefficient j3 has a maximum 
value. For even larger value of uj^, the conditions (25) 



and ( 26 1 are no longer fulfilled and both coefficients van- 
ish. It is energetically unfavorable to flip the spin of an 
electron against the direction of the effective constant 
magnetic field B^ s . As a consequence, as u>l increases, 



more energy is needed to flip the electronic spin to the 
direction opposite that of the field. This causes a to de- 
crease with increasing ojl- Additionally, the larger the 
ratio Wi/r, the less probable it is that spin-angular mo- 
mentum will be exchanged between the molecular spin 
and the itinerant electrons. For ujl = 0, the molecular 

spin is static, i.e. S = 0. In this case T(t) — 0. Ac- 
cording to Eq. (23 1, the coefficient a drops to zero for 
values lol = ±(e-f — ej,), i.e., for values \<jJl\ equal to the 
Zeeman splitting energy. At u>l = — ( e t — e i) = ~0-4eoi 
the level — wl is aligned with the Zeeman level e-j-, 
while the level c.^ + is aligned with the Zeeman level 
£[.. Therefore, there are only two channels available for 
electronic transport. At this point, the Gilbert damp- 
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FIG. 9: (Color online) (a) Gilbert damping coefficient a and (b) coefficient /3 as functions of the precession frequency ujl = 
g\isB c of the spin S of the SMM, with B c = B c e z , for different tunneling rates F at zero temperature. The applied bias voltage 
is eV = fix, — fiR = 2eo, with = 0. The other parameters are the same as in Fig. 7. In the case of the smallest F (red lines), 
the coefficient a is positive for e-|- < eV < + uil, ujl > , and negative for ej_ < eV < — uil, u>l < 0. The coefficient /3 has 
two local maxima and two local minima corresponding to the resonance of [il with energy of the levels in the QD. 



ing vanishes while the coefficient /3 reaches an extremum 
value for large values of T. With further decrease of uil, 
the relative position of the spin-up levels and — wl 
is reversed, as well as the relative position of the spin- 
down levels ei and e-j- + ojl- In this case, the transfer 
of spin angular momenta occurs in such a way that the 
coefficient a takes negative values and hence the molec- 
ular spin S becomes unstable. For lo l < and T <C \ui L \ 
(red lines), at the value of ui^ for which fi L = — lu^, 
the coefficient a has a negative minimum value while the 
coefficient /3 has a negative local minimum. The coeffi- 
cient a then increases with a further decrease of u>l as 
long as < /i^ < — lj^ ■ At the value of for which 
fjtjj = ej_, a drops to zero while (3 has a small local max- 



imum. According to Eq. ( 23 1 , the Gilbert damping also 
does not occur for wj, = — ei, which is realized if S is 
perpendicular to B c . 



IV. CONCLUSIONS 

In this paper, we have first theoretically studied time- 
dependent charge and spin transport through a small 
junction consisting of a single-level quantum dot coupled 
to two non-interacting metallic leads in the presence of 
a time-dependent magnetic field. We used the Keldysh 
non-equilibrium Green's functions method to derive the 
charge and spin currents in linear order with respect to 
the time-dependent component of the magnetic field with 
a characteristic frequency u>. We then focused on the case 
of a single electronic level coupled via exchange interac- 
tion to an effective magnetic field created by the preces- 
sional motion of an SMM's spin in a constant magnetic 



field. The inelastic tunneling processes that contribute 
to the spin currents produce an STT that acts on the 
molecular spin. The STT consists of a Gilbert damping 
component, a, as well as a component /?, that acts as an 
additional effective constant magnetic field and changes 
the precession frequency lul of the molecular spin. Both 
a and f3 depend on uj, and show a non-monotonic de- 
pendence on the tunneling rates, T. In the weak cou- 
pling limit r <^ Wi, a can be switched on and off as 
a function of bias and gate voltages. The coefficient (3 
correspondingly has a local extremum. For T — > 0, both 
a and f3 vanish. Taking into account that spin trans- 
port can be controlled by the bias and gate voltages, 
as well as by external magnetic fields, our results might 
be useful in spintronic applications using SMMs. Be- 
sides a spin-polarized STM, it may be possible to detect 
and manipulate the spin state of an SMM in a ferro- 
magnetic resonance experimentPsHU] a nd thus extract in- 
formation about the effects of the current-induced STT 
on the SMM. Our study could be complemented with 
a quantum description of an SMM in a single-molecule 
magnet junction and its coherent properties as these ren- 
der the SMM suitable for quantum information storage. 
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